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ABSTRACT 

The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. 

Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and 
grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & 

Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a 
compact emission region in the outer disk indicating a strong concentration of big dust grains. Tracing smaller 
grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to 
the north-west of the central star, which matches the VLA clump position. The compactness of the Ka emission 
is expected in the context of dust trapping, as big grains are trapped more easily than smaller grains in gas 
pressure maxima. We develop a non-axisymmetric parametric model inspired by a steady state vortex solution 
with parameters adequately selected to reproduce the observations, including the spectral energy distribution. 
Einally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially 
coincides with a spiral-like feature seen in scattered light, while the VLA clump is offset from the scattered light 
maximum. Moreover, the ALMA map shows a decrement that matches a region devoid of scattered polarised 
emission. Continuum observations at a different wavelength are necessary to conclude if the VLA-ALMA 
difference is an opacity or a real dust segregation. 

Subject headings: protoplanetary disks - planetdisk interactions - techniques: interferometric 


1. INTRODUCTION 

The formation of terrestial planets or of giant planets by 
core accretion requires the build up of a population of plan¬ 
etesimals. To form these bodies, dust grains need to grow 
from micron to kilometer sizes, i.e. grains need to increase 
in size by at least nine orders of magnitude in size scale, un¬ 


(Safronov & Zvjagina[ 

19691 [Pollack et al.[ 

[1996[). This pro- 

cess is a race against t 

he disk dispersal process ([Hollenbach[ 


star ([Hartmann et al.|1998|). 

As dust grains grow, they start to decou ple from the gas 


and d rift radially towards the central star ( jWeidenschilling 
1911) . Therefore, planetesimal formation has to happen in 


time scales shorter than radial drift. One way to halt the in¬ 
ward drift is by developing a local maximum in the radial sur- 


in density maxima ([Whipple] 19721 [Barge & Sommeria|1995 ; 

Haghighipour & Boss[[2003l 

[Chian^ 

1 8l Youdin[[2010), espe- 

daily big grains ([Brauer et a] 

[.|2008t[Pinilla et al.|2012). Dust 


of gaps and cavitie s, commonly called transitional disks ( Es 
jpaillat et al.|2014| ). These gaps and cavities are also detected 
as infrared dips in the spectral energy distribution (SED). In 
these dust traps, and depending on the level of gas turbulence, 
grains can accumulate for periods comparable to the disk life¬ 
time. Thereby benefiting grain growth, and prod ucing an en- 
hancement in the population of larger dust grains (jBrauer et al. 
|2008l|PimII5 et al.|2012D . 

smarino @ das .uchile. cl 


dust trapping ( 

Barge & Sommeria[ 

119951 ITanga et al.||1996 

Johansen et al. 

2004 

llnaba & Barge|2006l Lyra et al.||2008 

2009[ Regaly et al.[ 

2012[ Lyra & Lin 20l3j). Large scale 


the tr ansition disks a round SAP 206426 , SR 21 ( jPerez et ^ 
[2014] ), LkHo f 330 (jlsella et al.||2013|), and particularly in 
HD 142527 (jCasassus et al.||2013Er Fukagawa et al.||2013| ) 
and IRS 48 ( [van der Marel et al.|2013 ), where a stronger than 


usual azimuthal contrast is observed. These observations sug 
gest dust trapping by large-scale vortices. In this scenario it is 
expected that the dust spatial distribution will depend strongly 
on the grain size, as aerodynamic drag depends on the grain 
cross section. Big grains will concentrate in more compact 


regions than small grains (Barge 

; & Sommeria 


Klahr 

& Henning|1997[[Dubrulle et al.l 


[Inaba & Barge 

2006[ 

Meheut et al.[[2012[ Lyra & Lin] 

2013[ Raettig et al. 2015J. 


dust continuum emission should be more compact than at 
shorter wavelengths, as it traces the mass distributio n of grains 


with sizes comparable to the observed wavelength (Beckwith 
jet al.||199^ |Testi et al.||2003l |Bimstiel et al.|[2QT3 ). On the 


other hand, [Mittal & Chiang| ( |2Q15| ) show that large and sim 
ilar structures to the observed asymmetries in protoplanetary 
disks can be formed by gravitational global modes, which are 
currently difficult to discriminate from vortices. Thus, in the 
absence of tracers of the gas mass, it is crucial to use multi¬ 
wavelength observations to diagnose if dust trapping around 
a pressure maximum is present. 

In this paper we study the transition disk around the Herbig 
Ae star MWC 758 (HD 36112). The main parameters for this 
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syst em are given in Tab le [T] The central star is 3.7 ± 2.0 Myr 
old ( |Meeus et al.|2012| ), and optical spectroscopy reveals vari¬ 
able pronies'TuggesnngJet-likemtiomogene^s at different 
distances from the star ( [Beskrovnaya et al. 1999 ). The revised 
Hipparcos parallax d ata put the star at a distance of 279^53 pc 
( van Leeuwen|2Q07) ), although most previous st udies use the 


old Hipparcos estimated distance of 200_4o r 


cker et al.||1998| ). Studying the ^^CO(3-2) kinematics, 
determined a stellar mass of 2.0 ± O. 2 M 0 


Isella 
and a 


et al. (201 

disk orientation with an inclination (i) of 21 ± 2° and a po 
sition angle (PA) of 65 ± 7°, all consistent with Keplerian 
rotation. Despite its bright infrared excess, a cavity of 0.37" 
(73 astronomical units (au) at 200 pc) has been inferred from 
dust millimeter emission suggesting a depletion of millimeter 
grains in the inner regions of the disk possibly due to the pres¬ 
ence of a low-mass c ompanion within a radius of 42 au (using 
the revised distance) ( [Isella et al.|2010tpkndrews et al.|2011| ). 
In addition, the S EP shows a 10 /xm dip consistent with a pre- 
transitional disk ( [Grady et al.||2013| ), but with a strong near- 
infrared emission coming from sub-au scales observed with 
VLTI ( [Isella et al.|20Q8l ). 


TABLE 1 

Main parameters oe MWC 758. 


Parameter 

Value 

Reference 

RA (J2000) 

05^ 30"^ 27.530^ 

(1) 

DEC (J2000) 

+25° 19' 57.082" 

(1) 

Age 

3.7 ±2.0 Myr 

(2) 

Stellar Mass 

2.0±0.2Mo 

(3) 

V 

^21 

(4) 

B-V 

0.3 

(4) 

H 

6.56 

(5) 

Distance 

279+lt pc 

(1) 

Disk inclination 

21 ± 2° 

(3) 

Disk PA 

65 ±7° 

(3) 


^ (l)lvan Leeuwen|20071 (2)|Meeus et al.|20121 (3)|Isella et al.|2010|; (4) 
[ H0g et al. [ 1000[ (5) [ cjutrl et al. j lOO^ -"-' '-"-' 


The disk structure shows deviations from an axisymmetric 
disk. Spiral arms are detected from near-IR scattered light and 
thermal emission images extending from the sub-millimeter 
cavity to the outer disk ( [Grady et al. 2013t [Benisty et al.[ 
201 5[), possibly p roduced by a massive planet at ^160 au 
)ong et aL] ( [2015[ ). The drop in the polarised intensity beyond 
the spiral arms can be interpreted as a shadowing effect by the 
arms, as the disk extends farther out in the millimeter. In ad¬ 
dition, at millimeter wavelengths a peak intensity in the dust 
continuum has been detected to the north-west of the central 
star after subtraction of a best-fit axisymmetric model, sug¬ 
gesting an asymmetry in the millimeter dust gra in distribution 
in the outer parts of the disk ( [Isella et al.|2010 ). Moreover, in 
the same study, asymmetries m CO emission were observed 


that could b e due to a warped optically thick inn er disk (Isella 
201Q[ ), similar to the case of HD 142527 ( [Marino et al. 


et al. 


2015 ). If the southern spiral arm is trailing, the disk is rotating 
with a clockwise direction in the plane of the sky, where the 
north side of the disk is the nearest. 

We analyse new VLA Ka & Ku observations and we com¬ 
pare with ALMA archival data to study the distribution of dust 
grains and search for evidence of grain size segregation due 
to dust trapping. In Section]^ we describe the data and imag¬ 
ing. In Section we investigate the origin of the VLA Ka 
& Ku asymmetric emission, comparing with archival ALMA 


and SPHERE data. In Section 3.4 we describe a paramet¬ 
ric asymmetric disk model inspired by a modified vortex so¬ 
lution, that can reproduce the basic morphology seen in the 
VLA and ALMA images and the SED, to try to interpret the 
observed asymmetries. In Sectionwe discuss about the ori¬ 
gin of the asymmetries: A large seme pressure maximum, e.g. 
a vortex formed at the outer edge of the cavity producing dust 
trapping, or an embedded massive forming planet accreting 
material while heating its its surroundings. Linally in Section 
Owe summarise the main conclusions of this work. 


2. OBSERVATIONS AND IMAGING 

In this section we describe the observations and imaging of 
the VLA Ka (-33 GHz) & Ku (-15 GHz) and ALMA Band 
7 (—337 GHz) data. The image synthesis was carried out us¬ 
ing a non-parametric least-squares modeling technique with 
a regularization term called entropy from the family of max¬ 
imum entropy methods (MEM). We call the whole algorithm 
and code uvmem and the resulting deconvolved model images 


omy can be found in Pantin & Starck[ 

([1996); Casassus et al. 

([2006[);[Levanda & Leshem ([2010[);[C 

"asassus et al. ( 

20l3bf; 

Warmuth & Mann (2013 [); [Goughian & Gabuzda[([20] 

3b. It is 


beam derived from the response of the algorithm to a point 
source with the same noise level as the observations. The 
MEM synthetic beam is usually finer than the standard Clean 
beam with natural or hriggs weighting. 

The deconvolved model images are ’restored’ by adding the 
dirty map of the residuals of uvmem and convolving with 
a Clean beam characteristic of natural or briggs weighting. 
These images can be compared with Clean images and are 
labelled as ’restored images’. 


2.1. VLA 

The VLA observations form part of the project VLA 13B- 
273. The MWC 758 observations were acquired on 6 runs: 3 
nights in October and November of 2013 in Ka band and an¬ 
other 3 nights in January of 2014 in Ku band. In Ka band the 
target was observed in B 0 configuration with minimum and 
maximum projected baselines of 190 m and 10.5 km, while 
in Ku band it was observed in BnA 0 configuration with min¬ 
imum and maximum projected baselines of 80 m and 23 km. 
The array included 27 antennas of 25 m diameter, and the total 
integration time on source was Ih 21 min in each band. 

Lor Ka band, the VLA correlator was configured to produce 
64 spectral windows from 28.976 GHz to 37.024 GHz, with a 
spectral window bandwidth of 128 MHz, each divided into 64 
channels to study the continuum emission with a resolution 
of 2 MHz. In Ku band, the correlator was configured simi¬ 
larly with 64 spectral windows covering from 11.976 GHz to 
18.224 GHz, each divided into 64 channels of 2 MHz width. 

In all the observations, 3c 138 was used as primary fiux cal¬ 
ibrator, while 10547-^2721 was used as phase calibrator with 
4 observations of Imin 30s between the target observations of 
9 min. After looking at the resulting images from each run 
to check the quality of the data, we excluded one of the three 
observing runs with Ka band. 

2.1.1. VLA images 

In Ligurel^we present the restored VLA Ka image and re¬ 
stored image after subtracting the central component with a 
point source fit to the visibilities. The Ka image presents 
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Fig. 1.— Restored VLA Ka 33 GHz) images. Left panel: restored image. Right panel: restored image after subtracting the stellar emission with a point 
source model to the visibilities. The beam size in the Ka images is 0.23" x 0.22" and is represented by a white ellipse in both images. The x &y axes indicate 
the offset from the stellar position in RA and DEC in arcsec, i.e. north is up and west is right. The stellar position is marked with a red star. 


a noise level (a) of 6.6 /iJy beani“^ and the beam size is 
0.23" X 0.22" with natural weights. Two bright compact re¬ 
gions separated by 0.5" stand out from some extended emis¬ 
sion, which is not very well resolved due to low signal-to- 
noise. The southern compact emission has a peak intensity 
of 63.0 /iJy beam“^ and is located at the position where the 
star is e xpected according to th e stellar coordinates and proper 
motion ( [Perryman et al. 1997) . The other compact emission is 
44.0 / iJy beam~^ and is located to the north-west of the star, 
where jlsella et al.|(|2010|) found an intensity maximum at mil- 
limeter wavelengths. To the south of the star there is also disk 
emission with a peak intensity of 30.0 /iJy beam“^, slightly 
higher than 4cr. This probably correspond to dust thermal 
emission from the disk and puts an upper limit to the disk 
emission at this frequency. However, the morphology of this 
emission is hard to elucidate given the noise levels and reso¬ 
lution. The northern clump is separated from the star by 0.5", 
which at a distance of ^280 pc translates to a de-projected 
radius of ^150 au in the disk plane. The total flux at 33 GHz 
inside a circumference of 1.0 arcsec radius is 366±49/iJy and 
289 ± 33/iJy after subtracting the stellar emission. 

The Ku data only shows compact emission coming from 
the star location with a total flux of 67.1 ±7.7 /iJy and an im¬ 
age noise level of 4.1 /iJy beam“^. The synthetic beam of 
this observation corresponding to natural weights has a size 
of 0.56" X 0.18". We are not showing the Ku image because 
it is featureless and only shows a point source at the stellar 
position. 


2.1.2. Flux loss 

We neglect flux loss in the VLA data. To study this possibil¬ 
ity in the VLA Ka observations due to missing baseline spac¬ 
ing and limited u-v coverage, we simulate observations using 
the VLA Ka u-v coverage and a model im age o f the disk at 
337 GHz (the model is described in Section [T4| ). At this fre¬ 
quency the disk emission is more extended than at 33 GHz, 
and thus flux loss is more likely. We scaled the intensity lev¬ 
els of the model image at 337 GHz to have a simulated peak 
intensity equal to the Ka peak intensity. Gaussian noise was 
added to the model visibilities to reproduce the same noise 


level than the restored VLA Ka image. The flnal simulated 
image has a total flux of 620±27 /iJy inside a circumference 
of 1.0 arcsec radius, centered at the stellar position. This is 
consistent with the scaled model flux of 609 /iJy within the 
estimated errors. Therefore, we conclude that there is no flux 
loss due to missing uv sampling in the case of the VLA data. 

2.2. ALMA 

To study the disk continuum emission at millimeter wave¬ 
lengths we made use of ALMA archival data in Band 7 of 
MWC 758 from the project 2011.0.00320.S ( |Chapillon|20T5] ). 
The observations were conducted over two mgEfsTiTAu^t 
and October of 2012. The ALMA correlator was set in Fre¬ 
quency Division Mode (FDM) to provide 4 spectral windows 
divided in 384 channels, centered at 337.773 GHz with a to¬ 
tal bandwidth of 234.363 MHz, 344.469 GHz with a total 
bandwidth of 234.363 MHz, 332.488 GHZ with a total band¬ 
width of 233.143 MHz and at 330.565 GHz with a total band¬ 
width of 233.143 MHz. The observations were carried with 
29 antennas of 12 meter diameter. The minimum and max¬ 
imum projected baselines were 15 and 375 m and the total 
time on source was 12 min. We flagged the lines CO J =3- 
2 (337.79599 GHz) and J =3-2 (330.58797 GHz) to 
study the dust continuum emission. 

In Figure we present the ALMA MEM image using 
uvmem and a restored image that can be compared to a Clean 
image. The model image has an approximate resolution of 
0.31" X 0.18", derived from the response to a point source. 
The restored image was produced by convolving with a Clean 
beam of size 0.64" x 0.40" characteristic of hriggs weighting 
to achieve a better resolution. In the restored image, the total 
flux is 210 ± 5 mJy, with a peak intensity of 55.0 mJy beam 
and a noise level of 1.3 mJy beam“^. 

In Table [2] we summarise the new flux values of MWC 758 
added by this work along with global spectral trends. 

3. ANALYSIS 

The VLA Ka and ALMA maps show that the disk of 
MWC 758 is far from an axisymmetric disk. The disk shows 
a peak of intensity to the north-west side of the disk, centered 
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Fig. 2.— ALMA maps at 337 GHz (Band 7). Left panel: MEM non-parametric model (regularized with the maximum entropy method) with an approximate 
angular resolution of 0.31" x 0.18". Right panel: Restored image adding residuals and convolving with a Clean beam corresponding to briggs weighting 
(0.64" X 0.40"). The respective beams are represented by blue ellipses. The x & y axes indicate the offset from the stellar position in RA and DEC in arcsec, 
i.e. north is up and west is right. The stellar position is marked with a red star. 


TABLE 2 

VLA AND ALMA FLUXES 


Measurement 

Value 

ALMA band 7 flux (337 GHz) [mJy] 

205 ±5 

VLA Ka (33 GHz) flux [/xJy] 

366 ± 49 

VLA Ka disk flu?0[/iJy] 

289 ± 33 

23(j^337 GHz spectral index 

3.4 ±0.1 

337-33 GHz spectral index 

2.8 ±0.1 


^Flux inside a circumference of 1.0 arcsec radius after subtracting the cen¬ 
tral component _ 


Chapillon et al 


2008 


•^Data from ( 

at a distance of 0.5" fro m the star, consist ent with previous 
millimeter observations ( jlsella et aL]|2QlQ| ). In both ALMA 
maps and more evident in the MEM model, a local intensity 
minimum appears at the center of the disk, confirming the 
presence of a cavity of millimeter grains. Unlike the case of 
HD 142527 and IRS 48, the maximum and minimum of in- 
tensit y are not opposite in azimuth (azimuthal wavenumber m 
= 1) (|Casassus et al.||2013at |Fukagawa et al.||2013[ |van der 
|Marel et al.|2013| ).~ Moreover, a second maxinium m azimuth 
appears in the ALMA MEM model to the south-west of the 
star, which translates to an azimuthal extension to the south 
of the northern intensity maximum in the restored image. 


ilar unresolved emission has been detected with ATCA obser¬ 
vation at the location of the central star in HD 142527 with 


a slightly higher spectral index of L0±0.2 (Casassus et al. 


On the other hand, the northern compact emission has an 
inter-band spectral index of = 3.1 ± 0.1 when we com¬ 
pare with the intensity maximum of the restored ALMA im¬ 
age. Using the predicted flux of the northern clump at 
15 GHz is ^6.0 /iJy, slightly above the noise level of the Ku 
data. 

Given the unknown temperatures, the uncertainties in the 
clump flux, the non detection at 15 GHz, and the lack of spa¬ 
tial resolution in the ALMA data, it is not possible to con¬ 
clude on the spectral index /3 for an optical depth power law 


t{u) (X z/^, preventing us from studying a direct relation 
between the observations and a dust size distribution in the 
northern clump. Additional observations with similar resolu¬ 
tion levels are necessary to derive a P index to directly study 
the dust size distribution in the clump. However, in Section 
|3.4| we implement a simple parametric vortex model for the 
disk that reproduces the SED and morphology of the VLA 
and ALMA maps. We use this model to interpret the VLA 


clump in Section 3.4 


3.1. Spectral trends 

To investigate the origin of the emission in the VLA maps, it 
is useful to compute the intra-band and inter-band spectral in¬ 
dexes, defined as = log(/ 2 //i)/log(z/ 2 /z^i). The Ka-Ku 
inter-band spectral index (afl) of the emission coming from 
the stellar location is 0.36 ± 0.20, while the Ku intra band 
spectral index is 0^15 = 0.8 ± 0 . 6 , computed from a point 
source fit to the Ku data at the different channels. Both val¬ 
ues are consistent with free-free emission from the central star 
associated with a stellar wind or stellar accretion. Assuming 
a spectral index between 0.5 and 1.0 from free-free emission 
and the VLA fiux level at 33 GHz, we expect emission lev¬ 
els of 02-0.1 mJy beam“^ at 337 GHz, well below the disk 
emission and the noise level reached in the ALMA data. Sim¬ 


3.2. Comparison between ALMA and VLA Ka maps 

In Figure we present an overlay between the restored 
VLA Ka image and the ALMA band 7 MEM model. The 
VLA Ka northern peak intensity matches the location of the 
ALMA band 7 maximum. At 337 GHz we expect that most 
of the emission comes from ^millimeter-sized grains. If this 
intensity maximum is tracing a maximum in the dust density 
distribution of millimeter-sized grains, the VLA and ALMA 
maximum could be due to dust grains being trapped in a 
pressure maximum of the gas. In the dust trap scenario, we 
would expect to observe at centimeter wavelengths a higher 
contrast between the intensity maximum and the rest of the 
disk and a more compact emission, probing the distribution 
of centimeter-sized grains that are more easily trapped in gas 
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pressure maxima. However, the ALMA maps are not as well 
spatially resolved as the VLA Ka image due to differences in 
the u-v coverage. In Sec. 33 we compare both datasets at the 
same resolution level. 

Figure]^ also shows disk emission to the south of the stellar 
position m the VLA Ka map that matches a second peak in¬ 
tensity in the ALMA MEM model. This peak is closer in than 
the northern clump, at an angular distance of ^ 0.3". How¬ 
ever this local maximum disappears in the restored ALMA 
image due to the larger beam. 



Fig. 3.— ALMA-VLA overlay. In blue the ALMA band 7 MEM model 
in an arbitrary color scale, while in red the restored VLA Ka image contours 
set at arbitrary levels to highlight the morphology (the lowest contour level 
is 4.2 cr). The beam size of the restored VLA VLA Ka image and ALMA 
MEM model is represented by a red and blue ellipse, respectively. The x & 
y axes indicate the offset from the stellar position in RA and DEC in arcsec, 
i.e. north is up and west is right. 


3.3. ALMA and VLA Ka map at the same resolution level 

To bring both datasets to comparable resolutions, we con¬ 
volved the restored VLA Ka image with the ALMA Clean 
synthetic beam of size 0.64" x 0.40" after subtracting the stel¬ 
lar emission. We call this map the degraded VLA image. In 
Figure]^ we show the contour levels 0.6, 0.75, 0.85 and 0.92 
peak intensity of the degraded VLA and restored ALMA im¬ 
ages. The degraded VLA image presents a morphology sim¬ 
ilar to the ALMA MEM model with two peak intensities: to 
the north-west and south of the star (see Fig. [^. The north¬ 
ern clump in the degraded VLA image is still more compact 
than in the restored ALMA image, with a larger contrast. The 
solid angle inside the 0.85 intensity maximum contour in the 
degraded VLA image is 0.09 arcsec^, while it is 0.23 arcsec^ 
in the ALMA map. This result supports the dust trap interpre¬ 
tation, finding that the 33 GHz dust emission is more compact 
than at 337 GHz. 

The peak intensity to the south of the star in the restored 
VLA Ka image remains after star subtraction and convolu¬ 
tion with the ALMA beam, with a morphology similar to that 
observed in the ALMA MEM model. The nature of this com¬ 
pact emission is not clear as it is just ^4 times the noise level 
in the restored VLA Ka image. At the ALMA resolution the 



Fig. 4.— ALMA-VLA map contours overlay. Blue contours: ALMA 
Band 7 restored image. Red contours: degraded VLA image (restored image 
after subtracting the star and convolving with the ALMA beam). The contour 
levels are 0.6, 0.75, 0.85 and 0.92 times the peak intensity of each map. The 
lowest contour level of the VLA map represents at this resolution. The 
dashed red and blue lines represent the 0.85 peak intensity level. The beam 
size of both maps is represented by a black ellipse. The x & y axes indicate 
the offset from the stellar position in RA and DEC in arcsec, i.e. north is up 
and west is right. 


cavity is hard to distinguish, although the peak intensity is off¬ 
set from the stellar position in both maps, suggesting a cavity 
depleted of big grains. 

To confirm the dust trap interpretation, ALMA observations 
in extended configuration are required to resolve the disk with 
a similar or finer resolution than the VLA Ka observations. 
Deeper VLA observations are also necessary to detect the rest 
of the disk. 


3.4. Parametric non-axisymmetric model 

We develop a parametric non-axisymmetric disk model, in¬ 
spired by the steady stat e vortex solution to the gas and dust 
distribution described in Lyra & Lin ( 201^ (hereafter L^13). 


The aim is to reproduce part of the morphology seen in the 
ALMA and VLA maps, and the global SED. The model con¬ 
sists of a central star surrounded by a disk of gas and dust. 


& Kuruczl 


radius of t 

Lb Rr. 


with a temperature of 8250 K and a stellar 
, to fit the SED in the optical with a visual 
1 [T998] ). We define 


extinction of 0.22 ( |van den Ancker etal 
two distinct disk zones (see Table |^: 


1. A dusty inner disk with small grain sizes ranging from 
0.1 to 1 /im, composed of amorphous carbon and sili¬ 
cate grains to reproduce the NIR excess. 

2. A non-axisymmetric outer disk composed of amor¬ 
phous carbon and silicate grains to reproduce part of 
the disk morphology seen in the VLA and ALMA im¬ 
ages, and the photometry from the SED in the FIR and 
radio. 


The total gas mass in the disk is 0.1 M^. The outer disk 
incorporates an analytic prescription for a steady state vortex 
defined for both gas and dust distribution, based on the work 
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TABLE 3 
Disk parameters 


Parameter 

Inner disk 

Outer Disk 

Tmin 

0.2 

20.0 

Tmax [SU] 

80.0 

300.0 

S(r) 

/ \ -0.5 

Eq.|^ 

H(r) [au] 

,, 

( 100.0 au) 


5.0 X 10-® 

7.1 XlO^ 

Mass fraction in Silicates 

0.9 

0.48 

^min Vp\ 

0.1 

1.0 

(^max [p] 

1.0 

4x10^ 


^The dust size distribution is n{a) oc a The dust internal density for 
silicate and amorphous carbon grains is 4.0 and 2.0 g/cm^ respectively. 


of L^13. We first define an axisymmetric gas bac kground in 


spired by the surface density distribution used in Isella et al. 
([20T0|, and following t he solution for the surface den sity ot a 
VISCOUS Keplerian disk ( |Lynden-Bell & Pringle|1974' ) 


U{r) = Ee ( - 


exp 


1 


2(2 - 7 ) 




2-7 

- 1 - 1 

r, 


-^(5+1) 


,(i) 

,( 2 ) 


where H is the scale height of the disk and S is defined as the 
ratio between the the Stokes number (St) and dt, a turbulent 
diffusion parameter that measures the strength of the turbu¬ 
lence in the disk. We assume a global 5t in the disk for sim¬ 
plicity, although the turbulence in the core of the vortex can 
be driven by different mechanisms than in the rest of the disk. 
We set Vc = 120.0 au and 7 = —2.0 to form an outer disk that 
matches the peak intensity and cavity in the ALMA observa¬ 
tions, and the SEP at long w avelengths, similar to the profile 
derived in Isella et al. ( |2010| ). However, these parameters are 
not well constrained due to the lack of spatial resolution in the 
ALMA data and SED degeneracies. 

We add the vortex steady state solution described by eq. 64 
in L^ 13 with a slight modification to include a global disk gas 
density field. We impose that the gas density described by the 
vortex solution has to match the gas background at the vortex 
boundary where the vortex flow becomes supersonic with re¬ 
spect to the disk. This happens when a (defined in L^ 13 as the 
axial elliptical coordinate corresponding to the vortex’s semi¬ 
minor axis) is equal to , where x is the vortex aspect ratio 

and ujy = /Vt is the dimensionless vortex frequency. The 
same applies to the dust, and a general equation for both gas 
and dust can be written as 


p(r,0,z)= p6(r,z)max{l, cpv-(r,0)}, 
pv = exp 


(r-ro)^^ roV 


2m 




= exp 


fix) 


2 x^^^ 


(^+ 1 ) 


(3) 

(4) 

(5) 


where c is the contrast between the density maximum inside 
the vortex and the background density field at the same 
radius (see Section]^ for more details), c is defined to have a 
global density field that matches the background density field 
at the vortex boundary. Hy = H / (fix) ^/Sf-l) is the vortex 
scale length, where fix) is a scale function defined in L^13 
and is a function of x- A smooth solution for the density field. 


considering a smooth shock perimeter can be written as 

p{r,(l),z,S) = pb(r,z)[l + {c-l)pvir,(l),z,S)]. ( 6 ) 

We set X = 5 (which is not well constrain, but between 
the range of accepted valu es to remove the elliptical i nsta- 
bility at the vortex centre, Lesur & Papaloiz^|2009| ) and 
6t = 3x10“^. To com pute ujy we use the G oodman- 
Narayan-Goldreich solution ( [Goodman et al.|1987] ). 

Using Eq. |^and reasonable values for dt = 10“^ — 10“^, 
it is impossib^ to produce a vortex with the contrast required 
to reproduce the asymmetry seen in the ALMA maps. To in¬ 
crease the contrast with the disk background, we artificially 
extend the vortex steady state solution beyond its original 
boundary to obtain a better match with the observations (see 
Sec. I^for a discussion on this). This translates into changing 


c to exp 




, where is the ratio between the 


increased vortex size and the vortex original boundary (see 
Appendix for a detailed description). We set = 3 to have a 
good match between the restored ALMA image and the model 
image convolved with the synthetic ALMA beam (see Eigures 

l^andl^). 

We implement the model using 9 dust species in total: 
2 dust species for the inner disk, representing silicates and 
amorphous carbon grains; 7 dust species for the outer disk, 
one accounting for small amorphous carbon grains (from 1 to 
10 yum), while the rest represent silicate grains from 1 yum to 4 
mm with different spatial distributions set by eq. The opac¬ 
ities for each dust species in the model are computed using 


the ’Mie Theory’ code written by jBohren & Huffman (1983) 



tions of the model images with the same u-v coverage and 
noise levels of the VLA Ka and ALMA observations. In ad¬ 
dition, we present the model SED that matches roughly the 
observations. As expected from the vortex model, at long 
wavelengths the disk emission is concentrated in the pressure 
maximum as it traces the distributions of bigger grains which 
are highly trapped inside the vortex. 

The simulated ALMA image shows a peak of intensity sim¬ 
ilar to the restored ALMA image, although the morphology is 
not exactly the same. An interesting difference is the second 
intensity maximum seen in the ALMA MEM image, which 
appears as an extension to the south of the maximum in the 
restored image. This cannot be reproduced by a single and 
unresolved density maximum as our model shows. A second 
vortex would be needed to account for this emission. In Eig- 
ure 1 ^ we present the model predictions and simulated obser¬ 
vations of a model with two vortices located at the position of 
the VLA northern and southern clumps. The southern exten¬ 
sion of the ALMA peak intensity is naturally reproduced, and 
also a decrement appears at the north-east side of the disk due 
to the elongated beam (see figure in Sec. [^in Appendix). The 
VLA simulated observations show two unresolved clumps at 
the location of the VLA northern and southern clumps. The 


^ http://www.ita.uni-heidelberg.de/dullemond/software/radmc-3ci/ 
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Fig. 5.— Parametric model predictions, (a) 33 GHz model image, (b) 
337 GHz model image, (c) Simulated VLA Ka observation using the model 
image, (d) Simulated ALMA band 7 observation using the model image, (e) 
Dereddened SED of MWC 758 (blue dots) c ompared with the model (red 
line). Photometric data points were taken from|H0g et al.|(|2000);|Cutri et al.| 
J2003k|Elia et al.M2005);|Chapillon et al.|(12Uuk);|rsella et al. 1^201 U) and the 
new VLA Ka ana archival ALMA band / data presented m this w ork. We 
deredden the obser vations assuming a visual extinction Ay =0.22 ^van den| 
[Ancker et al.|1998) and an total-to-selective extinction ratio of Ry=3.1. Ihe 
errorbars represent 2a errors. In (a)-(d) the a: & y axes indicate the offset 
from the stellar position in RA and DEC in arcsec, i.e. north is up and west 
is right. 


disk parameters are the same to the ones described above. The 
second vortex is located at a radius of 90 au and it has a size 
ratio Tg of 3.12. 

If the northern clump emission seen in MWC 758 in band 7 
and Ka data is produced by a dust trap, it is a possible location 
for planetesimal and planetary core formation, as grains can 
grow faster and avoid the inward drift. From our model the 
dust temperature at the clump position is ^30 K, but in the 
observations the brightness temperature of the clump at 33 
GHz is 1.0 K. Thus it is highly unresolved or it is optically thin 
at this frequency. From our synthetic model, the inferred dust 
mass inside the vortex boundary is ^ 45 enough mass 
to form several planetesimals that could lead to the formation 
of a planetary core of a gas giant. However, grain porosity is 


-16 -8 0 8 16 24 32 

Intensity [/xmjy beam”'] 


0 6 12 18 24 30 36 42 48 

Intensity [mJy beam”'] 


Fig. 6.— Parametric model predictions with two vortices, (a) 33 GHz 
model image, (b) 337 GHz model image, (c) Simulated VLA Ka observation 
using the model image, (d) Simulated ALMA band 7 observation using the 
model image. In (a)-(d) the a: & y axes indicate the offset from the stellar 
position in RA and DEC in arcsec, i.e. north is up and west is right. 


not included in our model, thus this inferred mass is probably 
higher than the real value. 

3.5. Comparison with SPHERE PDI data 

We also compare the VLA and ALMA observations with 
publicly available VLT SPHERE/IRDIS PDI (polarimetric 
differential imaging) data ( [Benisty et al.|2015] ). The data was 
taken on December 5, 2014 during the science verification 
time of SPHERE. MWC 758 was observed in DPI (dual po¬ 
larimetric imaging) mode in F-band (1.04 pm) with the IRDIS 
sub-instrument of SPHERE, with a 185 mas diameter coron- 
agraphic mask and a Lyot apodizer. 70 exposures of 64 sec¬ 
onds each were taken, of which 48 (total integration time 51.2 
minutes) were used for the reduction presented in this paper. 
PDI has proven to be a powerful technique to suppress the 
stellar light, and to first order approximates the refiected light 
off the disk surface by only measur ing the polarized fraction 
of the light (lAvenhaus et al.|2014 ). The frames were dark- 
subtracted, flat-fielded and bad pixels were removed, then 
centered. The polarized flux was calculated and instrumental 
effects acco unted for in a way simi lar to that described in the 
Appendix of |Avenhaus et al.| ( [2Q14| ). In Figure 0 the resulting 
SPHERE imageTevealsTEeaisFoutside of the inner working 
angle set by the coronagraph, including a strong, two-armed 
spiral. 

Unlike the case of HD 142527 , where a large cavit y of 
^130 au is see n in scattered light (Avenhaus et al.|2014|) and 
in the sub-mm ( [Casassus et al.|2Q13at|Fukagawa et al.|2013| ), 
exposing the inner rim of the outer disk to stellar radiation, the 
SPHERE PDI F-band image of MWC 758 shows that, down 
to the inner working angle (26 au), the cavity is not depleted 
of /im-sized dust grains traced at this wavelength. Moreover, 
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RA offset [arcsec] 

Fig. 7.— SPHERE/IRDIS-ALMA-VLA overlay. In grey scale, the po¬ 
larised intensity scaled by to highlight the surface density of the disk, 
obtained with SPHERE/IRDIS PDI. The decrement in the center correspond 
to the position of the coronagraph. The ALMA band 7 MEM map is repre¬ 
sented in blue contours, while the restored VLA Ka map is presented in red 
contours. The contour levels are arbitrary to emphasise the disk morphology. 
The X & y axes indicate the offset from the stellar position in the RA and 
DEC in arcsec, i.e. north is up and west is right. 


the polarized surface brightness drops abruptly be hind the spi¬ 
ral arm s, as if the arms shadowed the outer disk ( [Grady et~aL 
|2013|). The presence of fim size dust grains could be linked 
to flari ng of the gas disk, th e mass of a possible gap-opening 
planet ( jFouchet et al.|201()| ), or we might be looking at differ¬ 
ent stages of similar systems. 

In the same figure, contours from the ALMA band 7 MEM 
and restored VLA Ka image are overplotted in blue and red, 
respectively. No evident correlation appears between the most 
southern spiral arm and the dust emission at 337 and 33 GHz, 
although we lack spatial resolution in the ALMA data to trace 
fine structures. On the other hand, the disk emission to the 
south of the star in the ALMA and VLA maps matches part of 
the western arm, while the VLA northern clump maximum is 
slightly offset 0.2") from the SPHERE PDI maximum at 
the same PA. Such a radial displacement could be explained if 
there were a large wall of dust with a higher scale height shad¬ 
owing the disk at the clump radius. This is probably related 
to a moderate or low fiaring in the outer disk. 

The SPHERE image also shows a region devoid of scat¬ 
tered light emission to the north-east between two bright arms, 
matching also a decrement in the ALMA MEM model that ex¬ 
tends from the outer disk to the stellar position. The origin of 
this decrement could be related to a shadowing effect rather 
than a lack of material. 


could be linked to the vortex tail which also causes dus t trap¬ 
ping and can extend farther out from the vortex core ( Lesurj 


& Papaloizou|2009||2010[|Railton & Papaloizou||2014). Re 

Tfr&Ba 


cently, jSurvilfe & Barge ( |2015j ) showed that vortices larger 

than the scale height might exist. Moreover, numerical MHD 
simulations have shown that high contrasts in the gas density 
field can be reached when large vortices form at the edge of 
planet induced gaps, surviving for ~1000 orbits (jZhu & Stonej 

[ 2 ^ . - 

On the other hand, the azimuthal dust density structure in 
the outer disk is probably more complex than an azimuthal 
wavenumber m = 1. The second intensity maximum in the 
ALMA MEM model could be related to a second vortex at 
the inner edge of the outer disk. MHD simulations predict 
multiple vortices to form by the Rossby Wave Instability a t 


the inner edges of “dead zones” Lyra & Mac Low 
Moreover, the decrement to the north-east in the ALMA and 
SPHERE images cannot be easily explained by a decrement 
in the dust density and a single density maximum. 

An alternative explanation to the asymmetry in the dust 
continuum emission at 337 GHz is that there is a strong asym¬ 
metry in the dust temperature field in the outer disk, similar 
to HD 142527 where a large crescent shows two peak intensi¬ 
ties separated by a decrement i n temperature caused by shad¬ 
ows cast by a tilted inner disk ( jCasassus et al.|20f5| [Marino 
et al.[|2015[ ). If the possible decrement in MWC 758 is pro- 
duced by a tilted or warped inner disk we would expect to 
see clear decrements almost opposite in azimuth, similar to 


the decrements i n HD 142527 seen in polarized light (Aven 
haus et alT][2Q14[ ). An other possibility is that the spira 


seen m scattered light ( [Grady et al J20 13 [ [B enisty etal 


arms 


2015) 


cast shadows in the outer disk at different azimuth and radii. 
Moreover, a combination of both scenarios could be present. 
This c ould explain the asymmetric CO emission ( [Isella et ^ 
[2010[ ) as a temperature effect, the formation of spiral arms 
due to tidal forces from an inclined inner disk, or produced as 
a dynami c consequence of a perturbation of the gas tempera¬ 
ture field ( [Montesinos|2015[ ). However, the compact emission 
to the north-west of the star at centimeter wavelengths seen 
in the VLA maps suggests that the dust concentration is real, 
although azimuthal structure in the dust temperature field is 
expected as the disk is not axisymmetric. 

The compact emission in the VLA Ka data could also be 
explained by a companion object, e.g. a forming planet ac¬ 
creting material. This planet could be also the responsible for 
the observed spiral arms in scattered li ght, as it is close to the 
planet proposed by [Dong et al. ( [2015 ). Such forming planet 
would have higher temperatures than a passive disk as it is 
accreting, also heating the disk around it. ALMA multi-band 
continuum and gas line observations would be ideal to disen¬ 
tangle if the VLA and ALMA band 7 asymmetries are density 
or temperature effects and if the VLA clump is accreting or 
not. Additional H-a high contrast imaging could reveal the 
presence of an accreting companion. 


4. DISCUSSION 

MWC 758 shows spiral arms, non-axisymmetric dust ther¬ 
mal emission and asymmetries in the CO emission. The ori¬ 
gin of all the disk asymmetries is hard to elucidate. A pos¬ 
sible explanation for the non-axisymmetric disk seen at mil¬ 
limeter wavelength is a dust trap caused by a pressure max¬ 
imum. A steady state vortex as the one proposed in L^13 
cannot produce such a high contrast between the dust den¬ 
sity in the vortex core and the disk background. A solution 


5. CONCLUSIONS 

We report new VLA observations at 33 and 15 GHz obser¬ 
vations of the disk around the star MWC 758. The disk emis¬ 
sion at these wavelengths traces the distribution of centimeter¬ 
sized dust grains. We found a compact disk emission at the 
same location where previous SMA and CARMA observa¬ 
tions found d eviations from an a xisymmetric model at shorter 
wavelengths ( [Isella et al.|2010[ ). This compact emission is at 
0.5" 150 au) from the central star. 
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The VLA data also shows a strong compact emission com¬ 
ing from the star position, characterised by a spectral index of 
^0.5. This is likely due to free-free emission from the central 
star associated with a stellar wind or stellar accretion. 

We compare these data with ALMA archival data, produc¬ 
ing ALMA band 7 continuum deconvolved and restored im¬ 
ages with a finer resolution than previous SMA and CARMA 
observations. Both maps show non-axisymmetric disk emis¬ 
sion, with a peak intensity at the same location as the VLA 
clump, but more extended. These multiwavelength observa¬ 
tions fit in the context of dust trapping, as the bigger grains 
traced at lower frequencies are more concentrated around 
pressure maxima than the ALMA emission around the maxi¬ 
mum. However, in the ALMA maps the disk is not as resolved 
as in the VLA data. New ALMA long-baseline observations 
are needed to resolve better the disk morphology at millimeter 
wavelengths and to confirm the dust trapping effect. 

In addition, we compare these observations with a new 
reduction of archival VLT SPHERE PDI images. There is 
an offset of ^0.2" between the VLA clump and the north¬ 
western spiral feature. The same arm matches a second in¬ 
tensity maximum in the dust continuum seen in the VLA and 
ALMA maps. More interestingly, the decrement to the north¬ 
east in the ALMA MEM model matches a decrement in the 
polarized intensity, suggesting a shadowing effect rather than 
lack of material. 

We develop a parametric non-axisymmetric model for the 
disk, incorporating a steady state vortex solution based on 
L^ 13. Such a vortex model cannot explain the decrement seen 
to the north-east of the star and a second maximum seen in the 
VLA and ALMA data. To reproduce the peak intensity of the 
ALMA image, we required to extend the vortex solution artifi¬ 
cially beyond the shock perimeter. We also find that a second 
vortex could explain the extension to the south of the ALMA 
peak intensity and the decrement to the north-east as an image 
reconstruction artefact. 

MWC 758 is a very complex protoplanetary disk and there 
are still details to be addressed: what is the link between the 
possible vortex or vortices, and the spiral arms and apparent 
decrements? The missing piece of the puzzle might be found 
with a detailed study of the disk kinematics. 
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APPENDIX 


A. VORTEX MODEL 

Eollowing the procedure of L^13, to match the the vor¬ 
tex steady state solution from L^13 with a gas density back¬ 
ground, we impose both density profiles to match at the vor¬ 
tex boundary where the vortex flow becomes supersonic. This 
happens when a (defined in L^13 as the axial elliptical coor¬ 
dinate corresponding to the vortex’s semi-minor axis) is equal 
to where x is the vortex aspect ratio and ujy = Uy/U 
is the dimensionless vortex frequency. Assuming the vertical 
dependence is the same inside and outside the vortex, we can 
focus only in the midplane. The density field inside the vortex 
can be written as 


p{a, S) = cpb exp 


2iJ2 


(S + l) 


(Al) 


where c is a constant that has to be determined imposing the 
following condition 


pb=p{a = 


H 

X^v' 


,5), 


pb = cpb exp 


A¥(-5 + 1) 


(A2) 

(A3) 


finally obtaining c = exp + 1) 

solution is 


p(r, 0, z) = pb{r, z) max s 1, c exp 


Thus, the global 

a^fix) 


2iJ2 


-is + i) 


} 


A more realistic or smoother solution can be written as 

a^fix) 


p{r, (f), z) = pb{r, 2 :) 1 + (c - 1) exp 


2iJ2 


(A4) 

(5 + 1 ) 


(A5) 

If we extend the vortex from its original boundary to , 


c has to be redefined as c = exp 
final density field is described by Eq. [A5] 


VsH 

X^v ’ 

and the 


B. MEM VORTEX MODEL 

Eigure shows the MEM models of the ALMA simulated 
observations with a single and two vortices. 
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